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Abstract 

We have studied the effects of strong impurity scattering on disordered super- 
conductors beyond the low impurity concentration limit. By applying the full 
CPA to a superconductiong A-B binary alloy, we calculated the fluctuations 
of the local order parameters A^i, and charge densities, nA,nB for weak 
and strong on site disorder. We find that for narrow band alloy s-wave super- 
conductors the conditions for Anderson's theorem are satisfied in general only 
for the case of particle-hole symmetry. In this case it is satisfied regardless 
whether we are in the weak or strong scattering regimes. Interestingly, we 
find that strong scattering leads to band splitting and in this regime for any 
band filling we have a critical concentration where a superconductor-insulator 

quantum phase transition occurs at T = 0. 
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I. INTRODUCTION 



As is well known, s-wave superconductivity is possible even in highly disordered sys- 
tems. Examples are superconducting intermetalic alloys (such as AucSii^c, NbcSii^^, 
Moi-cRhM), heavily doped cubic perovskites (such as Bai^cKcBiOs with x ~ 0.4)1 and 
chevrel phases ( such as Cso.3MoS'2).@ What lies beyond these initially surprising facts is 
Anderson's theoremi according to which pairing of time-reversed states leads to a finite gap 
2A in the density of states. Namely, the one particle states involved in the pairing need 
not be eigenstates of any translation operators and hence both A and are only weakly 
influenced by the disorder. 

The key assumptions required for Anderson's theorem are (i) non- magnetic scattering 
only (time reversal symmetry) [ii) the self-consistent order parameter A does not fluctuate 
from configuration to configuration. Previously Gyorffy et a/.i examined the conditions 
for (ii) to hold, and found that spatial fluctuations in A could be neglected provided that 
coherence length 

tt/a 

is much greater than the lattice spacing, a. Anderson's theorem also follows from the classic 
theory of Abrikosov and Gorkov.l3i Their argument is based on perturbation theory and the 
proof requires that the real part of the self energy S(-E) varies slowly near the Fermi energy 
i^^p'.iS It turns out that this is equivalent to the condition ^ >> a in Eq. [^. 

The purpose of this paper is to examine disordered s-wave superconductors for which 
the spatial fluctuations of A cannot be neglected. The results of an such inquiry will be 
important for narrow band or short coherence length superconductors where A < Ep or 
^ ~ a, such as some of those mentioned above. In the work of Gyorffy et a/.i the order 
parameter fluctuations due to impurities were studied using a perturbative technique. Zhit- 
omirsky and Walkei0 also calculated corrections to beyond Anderson's theorem due to 
order parameter fluctuations evaluated perturbatively. In the present paper the goal is to 
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treat the spatially varying A and charge density fully self-consistently within the coherent 
potential approximation (CPA), and hence allow for the case of arbitrarily strong A or 
charge fluctuations. 

In this regime a number of interesting new issues arise. For example Ghosal et a/.lii 
showed that strong disorder scattering leads to strong spatial variations in A, with the 
formation of superconducting 'islands' where A is large and other regions where A is small. 
Moreover, they found that the spectral gap persists even when A is very small in large regions 
of the sample. By contrast it was argued by OppermanEl and Zieg lei0 that order parameter 
fluctuations lead to a finite density of states within the gap. Similarly, in the calculations 
of the non-self-consistent density of states by Annett and Goldenfeldlli an impurity band 
tail develops due to spatial fluctuations in A, and eventually this leads to gap-less s-wave 
superconductivity. Obviously the fluctuations in A can arise either due to randomness in 
the single particle site energy, at each atomic site, as in an alloy, or due to randomness 
in the attractive interaction potential, f/j, as considered by Litak and Gyorffy.El Here we 
consider only the former case of alloy type disorder. 

In this paper we calculate the spatially random, self-consistent order parameter and 
charge due to AB binary alloy type disorder. The on site energy is Sa on a fraction c 
of lattice sites and Eb on the fraction 1 — c. The Green functions are calculated using the 
coherent potential approximation(CPA)lll^. This has been shown to be exact in the limit of 
infinite dimensions@ and to reproduce the results of the self-consistent Born approximation 
for weak scattering {sa—Sb ^ W, where W is the bandwith ) and the self-consistent T-matrix 
approximation in the limit of low impurity concentration (c ^ 1). We calculate the self- 
consistent A or B site order parameter (A^, A^) and charge density (nA,^^). Interestingly, 
we find that particle- hole symmetry leads to A^ = A^, and we show that Anderson's 
theorem applies exactly in this case, even though ~ a or A ~ W. A similar result was 
obtained previously for weak (\U\ <^ and strong (|f/| ^ W)^ interactions, but only 

for the weak disorder limit. By contrast the result here is exact for both weak and strong 
scattering. Evidently, this latter result implies that ^ ^ a is not a necessary condition for 



Anderson's theorem. Furthermore in the extreme disorder hmit {ea — Sb > W) we have a 
disorder induced band-sphtting in the normal state. Remarkably, if the chemical potential 
lies in the band gap there is a superconductor to insulator quantum phase transition at 
T = 0. If the chemical potential lies inside one of the split bands, the normal state is 
metalic but there are impurity states inside the superconducting gap. 

Below, in section II, we describe the model, and our CPA formalism and report our 
numerical results for the cases of weak {ea — Sb « W) and intermediate {sa — Sb < W) 
scattering. In section III we show that the gap fluctuation vanishes in the case of particle- 
hole symmetry, whilst in section IV we argue that within CPA Anderson's theorem is exact 
in this case. Section V contain our results for the strong scattering {sa — Sb > W) case 
where we find a superconductor-insulator transition. 

II. THE MODEL AND CPA FOR S-WAVE SUPERCONDUCTING ALLOYS. 

We use an attractive U single band Hubbard model defined by the Hamiltonian: 



where cL, respectively, the creation and annihilation operators of electrons with 

spin a on the lattice site i, hi^ = c\^Cia- is the local occupation number operator, ^ibare is the 
chemical potential, and tij is the hopping integral from site i to site j. Ui is the attractive 
pairing interaction at site i. In all of the numerical calculations shown below the interaction 
potential is Ui = —3.2t. For other interaction strengths we found similar results, is the 
site diagonal random disorder potential which takes on values Ea with probablity c and Sb 
with probablity 1 — c. 

After applying the Hartree-Fock-Gorkov (HFG) approximation to Eq. ^jour task becomes 
a study of the Gorkov equation: 
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X G{l,j,lUn) = Sij 



(3) 



where G{i, j,tUn) is the Fourier transform, with respect to the complex-time variable r, of 
the Greens function G{i,j,T) = ^J2ujn^^'^"'^C^{i,j,iuJn).The self-consistency conditions for 
the local order parameter and charge density are 



P LOr, 



(4) 



and 



rii 



P Uln. 
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where is a positive infinitesimal. The task at hand is to solve the above equations for 
each configuration {si} of the site energy and pairing interaction, to obtain the local order 
parameter and charge {Aj, , Tij} and hence to calculate the configurationally averaged Greens 
function G{i,i;iuJn) = {G{i,i;iuJn',{Si}))- 

Whilst most of the salient features of disordered superconductors are well described by 
the Abrikosov and Gorkov(AG)@'i theory, recently the CPA has also been brought to bear 
on the problem.lill"0 Within the CPA one replaces the on site random potential with a site 
independent self-energy matrix 



V 



(6) 



which is determined by the condition that an A or B impurity, corresponding to Ea and Eb 
respectively does not scatter on the average. Thus the average Green function is given by 



E 
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X G^{l,j,lUn) 
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In the previous application of the CPA to the above mo the condition which 

determined S was implemented under the assumption that the pairing potential Aj does not 
fluctuate with the site energies. Therefore the self-consistency conditions (Eqs. ^ |^) were 
satisfied only on average. However, this presumption is not a necessary part of the CPA 
and, as will be seen presently, is unduely restrictive. Here, we consider the more general 
case where on an A or B site, with site energy ea or sb-, the pairing potential is allowed 
to be or Ab respectively and the two local gaps, A^ and A^, are determined by the 
condition that the corresponding local gap equations (Eq. ^ are separately satisfied. Thus 
if the probability that a site is occupied by an A atom is c and that for a B atom is 1 — c 
the generalized CPA condition is 

cfA{lUn) + (1 - c)fB{lLOn) = (8) 

where the single site T-matrices, Ta and Te, are given by 

fiiiu^) = V.iiun) (i - G'^il,j,tUn)V,iiUn)y' (9) 
in terms of the local, single site, scattering potential matrix 



Vi 



^ Si ~ Sii — fii — Aj — Si2 ^ 



(10) 



y —A* — S21 —Si — S22 + fJ-i J 
for i = A,B. 

Note that in the above expression not only Si and Aj are allowed to fiuctuate but the 
local /ij also take on different values on A and B sites. Evidently, such variations arise from 
charge fluctuations. In the present HFG approximation this is described by 

Clearly, once the above CPA problem has been solved for a set of Ea, Sb, ^a, ^Bt'^a and 
ub self-consistency requires that they are recalculated using the relations 

Aa,b = %^ E e'^-'^Gfi^{^, t; lu^) (12) 



nA,B = ^ E e^""''G^\•''(^, ^; ^o^„) (13) 

where G^'^{i,i]iuJn) is the Greens function matrix averaged over all configurations with an 
A or B atom, respectively, on the site i. In our calculation we assume, the non-magnetic 
case, that nj| = riii = ^rii, and hence 

fJ'A,B = fJ'bare — -UA,B'nA,B- (14) 

Moreover, the patially averaged Green functions G'^'^{i,i;iujn) are approximated by the 
Green function for A or B impurity in the CPA effective lattice described by the self-energy 
S. They are given by 

G^'^{i,i]tUn) = G''{i,i;tujn) + G''{i,i;tUn)TA,BG''{i,i]iuJn). (15) 

In short, equations 6 — 15 fully specify a self-consistent procedure, which when carried 
to convergence, constitutes the complete CPA for disordered superconductors, in the HFG 
approximation, for the model Hamiltonian in Eq. ^. Note that since the CPA is the mean 
field theory of disorder and the HFG approximation is that for superconductivity, the above 
theory should be regarded as the mean field theory which treats disorder and electron in- 
teraction simultaneously and on an equal footing. By treating the self-consistency only on 
average, earlier workjii'0@ did not include the effects of fluctuations in A and n. Previously 
the full CPA was only implemented to determine the influence of A. fluctuation on tJ In 
the remainder of this paper we investigate the full consequences of treating interaction and 
disorder together on the basis of Eqs. 6 — 15. 

Our results for the gap and charge fluctuations as functions of the average band filling 
n are shown in Figs.l and 2. To simplify matters we used a 2d square lattice with lattice 
constant a = 1 and band energy 

ek = —2t{coskx + cosky) (16) 

where t is the nearest neighbour hopping amplitude. In these calculations we have taken 
the energy difference 5 = — to be a significant fraction of the bandwidth W = 'it 



{S = 0.5t, 5 = 2t) and hence we are in the fairly strong scattering regime ( for 5 = 2t). As 
Fig.l shows 7^ except at the point n = 1, namely a half filled band. Clearly from 
Fig. 2 Ha ^ ub Sit any filling except or 2. Thus, unlike in previous calculations^ the A 
and n fluctuations are central features of our results. For emphasis we show in Fig. 3 the 
standard deviations of the order parameter 

= {{\SA,y\) = (|A?|) - (|A,|)2 = c(l - c){Aa - Asf (17) 

and charge density 

M„ = ((5n,)') = {nl) - {n,f ^ c(l - c){nA - nsf (18) 

as predicted by our CPA calculations. 

Remarkably, at half filling the fluctuations in the pairing potential go to zero, whilst the 
charge density fluctuations are at their strongest. To investigate the origin of this interesting 
phenomena we studied the case where c = 0.75 7^ 0.5 but the band is still half filled at fi = 1. 
For this case A^, A^ and A are shown in Fig. 4 as a function of temperature. Evidently for 
all T < Tc, A A 7^ A^ and hence 7^ 0. In what follows we unravel the root cause of this 
behavior. 



III. PARTICLE-HOLE SYMMETRY 

Recall that in Figs.l — 3 the order parameter fluctuations vanish, A^ = A^, for the case 
n = \ and c = 0.5, namely equal concentrations of A and B atoms. This special case is one 
where the Fermi energy is at the center of the band and the density of states is symmetric, 
and hence particle-hole symmetry occurs. For this case and only this case /i = 0, [ihare = 
and fiB = —fJ-A- Furthermore, since the self-energy always obey the same symmetries as the 
Green functions 

Sii(«u;„) = -S22(«cu„) 

Sl2(^^n) = T.l^{lUJn), (19) 
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for particle-hole symmetry there is also the property that 

Eii(zu;„) = -S*i(2a;„). (20) 



Namely, 



3?Sn(z^„)=0 (21) 



and consequently 



^Gl^iiiUn) = 0. (22) 



Noting that the CPA respects these symmetries we can rewrite the CPA condition in 
Eq. ^, as 

cV^' + (1 - c)V^' = G^i, I- lUn). (23) 



Combining this with particle-hole symmetry as described by Eqs. |19|, ^ |2^ implies that 



(A^ - Ab) + + (^A - ^) A^ + (/i^ - ^) A^^ = 0. (24) 

and hence, since the second bracket is non-zero, that 

A^ = As. (25) 
Therefore particle-hole symmetry implies the absence of the order parameter fluctuations. 

IV. PROOF OF ANDERSON'S THEOREM IN THE PARTICLE-HOLE 
SYMMETRIC DENSITY OF STATE CASE 

In non-magnetic disordered local s-wave superconductors the traditional argument lead- 
ing to Anderson's theorem assumes that the fluctuations of the order parameter are negligi- 
ble, A,; ^ A. Anderson's theorem shows that in this case there are no bound states between 
quasi-particles and impurity sites and therefore the quasi-particle energy gap is absolute. 
Namely, there are no impurity states inside of the gap.@ Furthermore is found simply by 
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replacing, the clean system density of states N{E), with its disordered system average N{E) 
in the gap equation. 

In an alternative route to the same result Abrikosov and Gorkovi'i use perturbation 

theory which implies that the self-energy can be approximated by 

i 

Sii(zu;„) = sign{ujn) (26) 

r 

where r is a wave vector and frequency independent quasi-particle life time. This is justified 
in the case of the non-self-consistent Born approximation by the assumption that the relevant 
energy scale is \uJn\ < {^d is Debeye frequncy) and hence only states near the Fermi 
surface are relevant. More implicitly it is assumed that near the Fermi surface the density 
of states is a constant and hence there is effectively particle-hole symmetry. 

However in general we cannot assume that Sii(2a;„) is of this form. In the case of narrow 
band superconductors, there is no Debeye cut off and so one cannot assume that only states 
near the Fermi level are significant. In particular 53Sii(zco'„) will not be a constant, and 
9ftSii(zci;„) need not be zero. Nevertheless, if for some reason particle-hole symmetry is 
obeyed, Anderson's theorem will obtain in full CPA, self-consistent Born, and T-matrix 
approximations. 

In Figs. 5 and 6 we illustrate using our explicit calculation the energy dependence of 
3?Sii(i?), QT,ii{E), ^{R{E)) and '^{R{E)) for the cases of particle-hole symetric or non- 
symmetric cases. Note that function R{E) is defined by 

R{E) = ^ [EniE + zO+) + Sn(-^ - ^0+)) , (27) 

and R{E) is the analytical continuation of 3?Sii(2a;„) to the real axis. Figs. 5 and 6 show 
that for particle-hole symmetry, n = 1 and c = 0.5, R{E) is equal to zero but in other cases 
R{E) 7^ 0. We shall now analyze the consequences of the particle-hole symmetry E —E. 

The Green function for CPA, self-consistent Born or T-matrix approximations can be 
written in the form 



k 



y A* ILOn + tk-flj 
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(28) 



where the renormahzed parameters A, u;„ and /i and rj^^ are given by 

Un = UJn{l ) 



/i = /i - 3?Sii(2a;„) 

,„ = (1 - (29) 

Note that these renormahzed parameters are the same as in the original paper of Abrikosov 
and Gorkovi, where S was computed in the non-self-consistent Born approximation. They 
found that SfJSn near the Fermi surface is independents of uj, and therefore eventually it can 
be absorbed into the chemical potential. The same procedure was followed Martin et al..^ 
However we proceed without these simplifications and retain the full energy dependence 
of ^T,ii{E) and $5Sii(E'). Now, particle-hole symmetry implies that 3?Sii(?ct;„) = 0, and 
therefore fi = fi. Using the above relations the gap equation becomes 

\U\ ^ p N{i, luj^) 



where 



Nii,tuJn) = ^Y.^Si~^~^~k) (31) 



and, efc = — is the renormahzed band energy. 



Surprisingly in the case of particle-hole symmetry the quantity A^(e, zoOn) in Eq. |3^ be- 
comes equal to the disorder average normal state density of states, namely, N{e) and hence, 
without further assumptions, Anderson's theorem obtains. To be quick and explicit we note 
that this last step follows from the property of a delta function that = a6{x), and hence 

N{~e, liOn) = ^ E -^(e - e~k) = = N{~e), (32) 

In short, we again have Anderson's theorem that Tc is given by the usual gap equation, but 
with the disorder average normal density of states. However, unlike other proofs, ours does 
not neglect the energy dependence of 3?Sii(i?) or QT,ii{E). 
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V. SUPERCONDUCTORS IN THE SPLIT BAND REGIME 



One of the main virtues of the CPA in the normal state is the fact that it describes 
band sphtting correctly.ill Namely, for ea — £b less than the bandwidth, 8t in our case, it 
predicts an effective band somewhere in between the bands of pure A or pure B metal, whilst 
for sa — bigger than the half bandwidth CPA predicts two, smeared, but nevertheless 
well defined, bands seperated by a gap. The two bands are at energies where there would 
have been an A or B band in one of the pure systems. In this later case the wave function 
corresponding to the A band is large mainly on the A sites and that in the B band is 
significant only on the B sites. In what follows we shall investigate superconductivity in this 
split band regime of CPA. 

To investigate the consequences for superconductivity in the above split band regime we 
have solved Eq.6-15 using a strong scattering potential {ea — ^b ~ W). From these solutions 
two interesting points emerged. Firstly we found that for the particle-hole symmetric case, 
reported in Fig. 8, = and hence even in this strong scattering state there are no fluc- 
tuations in A. Nevertheless, scattering has a large effect on the superconductivity through 
the configurationally average density of states N{E) in the Eq. |32|. Namely as shown in 
Fig. 8 A ^ for the critical strength of scattering 6 ~ 4.5t. 

The other interesting phenomenon is a superconductor to insulator transition. Recently 
Scalettar, Trivedi and Huscroftil discovered a superconductor-insulator transition in the 
disordered attractive U Hubbard model, using Monte Corlo simulations. They found that 
for strong disorder the superconducting gap is replaced by an insulating gap, both for weak 
and strong interaction U. Clearly the superconductor to insulator transition in our CPA 
calculaton has a similar origin, although our model of disorder is different. 

In general the band filling n is given by 

n = (1 — c)nB + cha- (33) 

where ua and ub are the partial averaged occupation numbers on A or B sites respectively. 
We shall discuss three different cases of Eq. case (i) (n^ = 0,nB < 2), case (ii) {ua = 
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0,nB = 2) and case {in) {ua 7^ 0,71^ = 2). In the first case the A band is empty and the 
filling of the B band(nB) is 

n , , 

nB = z . 34 

1 — c 

For the second case the B band is completely, doubly, occupied and the A band completly 
empty.In the third case the A band is partially occupied with 

n-2(l-c) 

riA = ^- -. (35) 

and the B band is fully occupied. The second case is a special case of Eq. Fig. 7 shows 
the A, B and average normal density of states for these three different concentrations. In 
Fig. 7(a) the A sites are approximately empty (n^ ~ 0) but there is more than one electron 
on the B sites. The graph plotted was for the case n = 1 and c = 0.25, and therefore from 
Eq. ^we have = |. Fig. 7(a) is a band metal with hopping of electrons from the B sites 
to the B sites. In Fig. 7(6) the A sites are almost completely empty, ua = 0, and the B sites 
are doubly occupied, = 2. Therefore there is a gap at the Fermi energy and the system 
is a band insulator. In Fig. 7(c) the B band is fully occupied while the A band is partially 
filled. For this case we had n = 1 and c = 0.75 and so by Eq. ^ the band filling of the A 
sites is Ua = §• Similarly Fig. 7(c) is a metal band with hopping from the A sites to the A 
sites. In Fig. 7(6) there is no hopping and we can regard this state as a band insulator. 

As illustrated in Fig. 8a, in the split band regime the superconductor gap closes, A ^ as 
6 — s> 4.5t but is replaced by an insulating gap for 6 > A.5t . Evidently, this can also happen 
in the non particle-hole symmetric shown in Fig. 86. The general condition for the 

superconductor-insulator transition is saturation of the B sites with 2 electrons of opposite 
spin while the A sites are empty. This happens when 

n = 2(l-c). (36) 

In this case the Fermi level lies in the gap and all the B states are filled and the A states 
are empty. For completeness we show the suppression of A by disorder for particle-hole 
symmetric and non particle-hole symmetric cases in Fig. 9, obeying the condition from Eq. |3^. 
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Now numerically we shall test our prediction, in Eqs. ^ |3^, for the two non-particle- 
hole symmetric cases: c = 0.25, n = 1 and c = 0.75, n = 1, as discussed at the beginning 
of this section. In this cases there is no superconductor-insulator transition. For a greater 
understanding of the details of the band splitting mechanism we plot the density of states, the 
order parameter and the charge density on A and B sites. In the first case the concentration 
of A sites is less than B sites, but in the second case the concentration of A sites is more 
than of B sites. As expected, in Fig. 10(a) we see that there are two split bands, a normal 
empty band on the A sites and a superconducting band on B sites. Conversely in Fig. 10(6) 
the B band is a normal doubly occupied band and the A band is the superconducting band. 

The effect of disorder on the average order parameter and the A and B sites local order 
parameters (A, A^, A^) and charge densities (n^, ub) are shown in Figs. 11 and Fig.l2. One 
can see that, for weak scattering the local order parameters and charge densities of A and 
B sites are approximatly the same but with increasing disorder the difference between them 
will increase, with one of A^ or A^ going to zero(depending on the relative concentrations). 
In Fig. 11 (a) and Fig. 12 (a) the concentration of A sites is less than that of B sites therefore 
A^ — i> 0, A^ — s> constant while ^ 0, ub ^ |. In contrast to this, in Fig. 11(6) and 
Fig. 12(6) the concentration of A sites is more than that of B sites, consequently A^ — >■ 
constant, Ab and ua — > |, ^ 2. 

Schematically Fig. 13 shows the physical mechanism of the band splitting in terms of 
A and B lattice sites. Clearly when all B sites become doubley occupied and all the A 
sites become empty there is no hopping and the system becomes insulating. Therefore as 
a function of the disorder, 6, or band filling n there is a T = superconductor-insulator 
quantum phase transition. This transition corroborates the suggestion of Scalettar, Trivedi 
and Huscroftii that disorder in superconductors can lead to a superconductor-insulator 
transition. However it shuld be noted that they used a uniform random distribution of site 
energies, Ei, wheras we used a binary alloy model. 



14 



VI. CONCLUSION 



In this paper, for the first time, we apphed the full CPA for disorder superconductors. 
In our calculations the self-consistency equations were solved fully within CPA so that self- 
consitency was properly satisfied on each type of CPA impurity sites, and not just on the 
average. It is in this form that the CPA-HFG is the only 'controled' mean field theory of 
disordered superconductors.0 

We plotted the local order parameter and charge density of A or B atoms of a binary alloy 
for both weak and strong scattering limits. We found that only for one special case are the 
order parameters of A and B sites equal, and consequently the fluctuations of A are zero. For 
this point we have shown analytically that the condition for Anderson's theorem is fuUfilled 
not only for weak scattering but also for strong scattering in a particle-hole symmetric band. 
By contrast neither, the density of states nor is exactly constant, although in the gap 
equation, the quasi-particle energy gap is absolute. In the gap equation the normal clean 
system density of states in the gap equation is replaced by the normal disorder average 
density of states, and therefore the only changes of the superconducting density of states 
and Tc come from the latter. 

In narrow band binary alloy s-wave superconductors, we showed that strong disor- 
der leads to two different interesting phenomena: (i) band splitting with a quantum 
superconductor-insulator phase transition at T = 0, (ii) band splitting without a phase 
transition. In this last case in terms of concentration and average band filling, one band 
is normal (doubly occupied or empty) and superconductivity is only present in the another 
band that is partially occupied. 
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FIGURES 




n 

FIG. 1. Self-consistent order parameters at A or B sites, A^, and the average A as a 
function of band filling. Here T = 0.008625t,c = 0.5 and S = 2t. Note that at half band filling 

A^ = Ab = A showing the fluctuation vanishes when there is particle-hole symmetry. 





n 

(I A? I) — lA'^p, and fluctuation of charge density, 



0.2 0.4 0.6 0.8 1 1.2 

FIG. 3. Fluctuation in A, Ma 
M„ = (|nf I) - as a function of band filling for c = 0.5 and T = 0.008625t. Note that 

the A fluctuation is zero at half band filling, while the charge density fluctuation is maximum. 
Therefore n = 1 obeys the conditions for Anderson's theorem: particle-hole symmetry and as a 
consequence of this there is absence of fluctuation in A. For other fillings this condition is not true. 




0.1 0.2 0.3 0.4 0.5 0.6 



T/t 

FIG. 4. Ayi, As and the average A as a function of temperature in the particle-hole asym- 
metric case for n = 1, c = 0.75 and 6 = 2t. Note that in this case A^, Ab and the average A go 
to zero at Tc ~ 0.525i. 
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FIG. 5. and Q'Eu{E) as a function of energy E. The case c = 0.5, n = 1 has a 

particle-hole symmetry and others are non-particle-hole symmetric. As it is obvious that for 
such a narrow band superconductors neither ^T,ii(E) nor QT,ii(E) are constant, unlike the weak 
scattering Born approximation limit. 
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FIG. 6. R{E) is the analytical continuation of ^Yiii(iojn) to the real axis. Here we show the 
real and imaginary part of R{E), ?R.{R{E)), Q{R{E)) for three different cases: the cases c = 0.5, 
n = 0.576, S = 2t and c = 0.25, n = 1, S = 2t are particle-hole asymmetric, and the case c = 0.5, 
n = 1, 6 = 2t is particle-hole symmetric. The temperature here is T = 0.008625i. Clearly only in 
the particle-hole symmetric case is R{E) = 0. 
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FIG. 7. Band splitting of the normal system at half band filling n = 1 and T = 0.008625t 
in the three cases: (a) c = 0.25, 5 = 8t. In this case conduction is in the B band, (b) c = 0.5, 
S = 10.lt in this case Fermi energy lies outside of both bands, therefore the state is an insulator, 
(c) c = 0.75, 6 = 8t. In this case the Fermi energy lies in the A band, so hopping occurs primarily 
in this band. 



20 




0.45 
0.4 
0.35 
W 0.3 
Zo.25 
0.2 
0.15 
0.1 
0.05 




1 


.1 1 1 

: 5=0.1t 


; (b) 


; \ \ 5=3.5t----- 




\ \ §=5t " 


1 





-15 -10 -5 5 10 15 

E/t 

FIG. 8. The superconductor-insulator phase transition in the strong scattering hmit for two 
different alloy concentrations and average band filling at T = 0.008625*. (a) A particle-hole 
symmetric density of states at n = 1, c = 0.5. (b) An asymmetric particle- hole density of states 
for c = 0.25, n = 1.5. Note that in both cases the superconducting gap closes and is replaced by 
an insulating gap for large 5. 
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FIG. 9. The fluctuation of A, Ma = (l^f I) — l^'^P; and the average A as a function of disorder 
strength, (a) The particle-hole symmetric case c = 0.5 and n = 1, (b) A particle-hole asymmetric 
band c = 0.25 and n = 1.5 at T = 0.008625t. Note that the A fluctuation in the particle-hole 
symmetric band is zero for all of disorder strengths. 
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FIG. 10. Density of states at A and B sites Na-, Nb and the average N at half band filling, 
n = 1 and r = 0.008625t, for the two cases: (a) c = 0.25, (5 = 9i and (b) c = 0.75, 5 = 9t. In 
the (a) case, A band is an empty normal band and superconductivity is only in the B band. In 
the case (b), the B band is a normal doubly occupied band and A band is a partialy occupied 
superconducting band. 
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FIG. 11. Order parameters A^, Ab and A at half band filling n = 1 in terms of disorder 
strength for the cases of: (a) c = 0.25 and (b) c = 0.75. For the case (a) and very strong scattering 
the order parameter on the A sites goes to zero. In case (b) the order parameter of the B sites 
goes to zero while that of the A sites goes to a constant. 
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FIG. 12. Charge densities ua and at half band fining n = 1 and T = 0.008625i in terms of 
the disorder strength for the case of: (a) c = 0.25, (b) c = 0.75. In the case (a) all of electons are 
in the B sites and A sites are empty in the strong disorder limit. The band filling of the B sites is 

= |. But in case (b) the B band becomes fully occupied while the A band remains partially 



occupied. 
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FIG. 13. The physical mechanism of the band splitting in the particle-hole symmetric case of 
c = 0.5, n = 1 for T = 0.008625t. The band splitting happens at approximately 4.5t. In this case 
all of the A sites become empty while the B sites become doubly occupied. Clearly this represents 
a band insulator. 
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